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Abstract. 

Due to the Fermi-Dirac statistics of electrons the temporal correlations of 
tunneling events in a double barrier setup are typically negative. Here, we 
investigate the shot noise behavior of a system of two capacitively coupled 
quantum dot states by means of a Master equation model. In an asymmetric 
setup positive correlations in the tunneling current can arise due to the bunching 
of tunneling events. The underlying mechanism will be discussed in detail in 
terms of the current-current correlation function and the frequency-dependent 
Fano factor. 



PACS numbers: 72.70+m,73.23.Hk,73.40.Gk,73.63.Kv,74.40.+k 
1. Introduction 

Resonant tunneling of electrons through bound states in mesoscopic double barrier 
devices still attracts wide research activities. Not only the current flow is of interest. 
Due to the granularity of the electron charge the tunneling current is a discrete 
stochastic process which leads to shot noise in the current signal pQ. The spectral 
power density gives information about temporal correlations which are not available 
by conductance measurements alone. In the pioneering theoretical work of L. Y. Chen 
et al. [2] a suppression of the zero-frequency spectral power density S with respect to 
the uncorrelated value Sp — 2el (Poissonian noise, e > is the elementary charge, I 
is the average current) in tunneling through double barriers was obtained. This refers 
to negative correlations in the tunneling current due to Pauli's exclusion principle: 
if the bound state becomes occupied by tunneling of an electron, then a subsequent 
electron "has to wait" until the state becomes empty. Hence, negative correlations are 
a fingerprint of the fermionic statistics of electrons. Contrarily, boson statistics lead 
to positive correlations, e.g. in the intensity-intensity correlation function of emission 
of photons from a black radiation source due to temporal bunching of emission events 
and references therein). 

For tunneling through quantum dots (QDs) Coulomb interaction generates an 
additional source of correlations. It was shown in Ref. that for an asymmetric 

setup of the tunneling rates in a system of two capacitively coupled QD states the 
current-voltage characteristic exhibits negative differential conductance. In this regime 
positive correlations in the tunneling current are present (super-Poissonian shot noise). 
Here, we give a detailed examination of this effect and show that its origin is related 
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to bunching of tunneling events due to Coulomb blocking. Positive correlations in 
tunneling current were also investigated in a double barrier resonant tunneling diode 
UJ, in coupled metallic QDs [7], and in a single-electron transistor setup |SJ- 



2. Model 



We consider two QD single-particle states with energies Exn which are coupled to 
emitter/collector contacts by tunneling rates ^ E / c (i — 1,2), respectively. The 
emitter/collector contact is assumed to be in local equilibrium and the occupation 
is given by the Fermi functions f E /c(E) and chemical potentials He/Ci respectively. 
An applied bias voltage eV = — fJ>c drives the QD system out of equilibrium. The 
tunneling energy E is given by the difference of energies before and after a tunneling 
process considering the charging energy U which is present if both QD states are 
occupied. The time evolution of the occupation probabilities P„ for the Fock states 
v = (ni,ri2) (n, S {0,1}, i — 1,2) is described by the Master equation P = MP. 
The 4x4- matrix M consists of the transition rates for tunneling into or out of the 
QD system which contain the occupation of the contacts and the tunneling rates (for 
details see 0). 

The stationary current is given by I = 2„ (j P ) with the stationary 

\=E / v 

occupation probability given by MP = and the operator for the current at the 
emitter barrier j , respectively. The current-current correlation function at the emitter 

barrier for t > is C EE (t) = (I E (t)I E (Q)) - I 2 = £„ (j T(t)j P°) + 5{t)el - I 2 

\— E — E / i/ 

with the time translation operator T(i) = expMi so that P(t) = T(i)P(0) holds 
|H| . Then, the spectral power density reads (Wiener-Khintchine theorem) See(u) = 
2 dtC ee (t) cos uit. As a measure for correlations at the emitter barrier the Fano 
factor is defined as a(u) = S ' E ^} I UJ ' > which is unity for uncorrelated tunneling, < 1 
for negative correlations, and > 1 for positive correlations in the tunneling current 
(super-Poissonian shot noise). Note that See(0) — Scc(0) = Sec(0) = Sce{0), but 
S(uj) ^ S EE (lo) for lu ^ 0. 



3. Results and Discussion 



We consider the case where both energies of the QD single-particle states are equal: 
Ei = P 2 i the tunneling rates to the emitter are equal: T E = T E = Te, and a 
bias voltage V is applied such that both single-particle states (1,0) and (0,1) can 
be occupied but the occupancy of the doubly-occupied state (1,1) is energetically 
forbidden. Then, the zero-frequency Fano factor is derived analytically to be 

, n s 1 . 2[ 7l 2 + 7 2 2 - 7i7 2 (2 + 7i+7 2 )] 

a(0) = l+ ■ ■ -2 (1) 

(71 + 72 + 7172) 

with the ratios of collector and emitter tunneling rates 7^ = T^/Ve {i = 1,2). The 
dependence of the Fano factor on 71 and 7 2 is depicted in Fig. la. The Fano factor 
is symmetric with respect to an interchange of 71 and 72. Two different regions 
corresponding to a < 1 and a > 1 are labeled by N and P, respectively. The region 
./V indicates negative correlations in the tunneling current for 7 j > f (i = 1,2). In 
this regime, mutual Coulomb blocking of the single-particle states (1,0) and (0,1) 
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Figure 1. a) Fano factor «(0) vs. 71 = F^/F^; and 72 = F^/Fe (white 
lines correspond to ct(0) = 1 and confine the regions P and N). b) Schematic 
realisations of current Ii vs. time t through level i = 1 (upper plot) and i = 2 
(lower plot) for 71 <C 72 =1. 



is negligible since the levels are mostly empty in time average due to the stronger 
collector coupling. If either one ratio or both ratios become smaller than unity the 
tunneling is positively correlated (region P in Fig. la) except for 71 ~ 72 where a 
approaches unity. Hence, the conditions for positive correlations are: 71 < 1 or 72 < 1, 
and 71 7^ 72- To illustrate the effect of positive correlations we consider a realisation 
of time-dependent currents Ji(i) and I%(t) through the levels 1 and 2, respectively, in 
Fig. lb. There, it is assumed that 71 <SC 72 =1- If an electron jumps from the emitter 
into level 1 (black peak) this level becomes occupied. Until the level becomes empty 
by tunneling of the electron out in the collector (grey peak) no electrons can jump in 
level 2 (Coulomb blocking). Only if level 1 is empty, a current I2 is flowing. Hence, a 
bunching of tunneling events in current ^(t) occurs. 

In Fig. 2 the frequency-dependent Fano factor a{oS) and in the insets of Fig. 2 
the current-current correlation function Cee if) at the emitter barrier are shown. The 
dashed curve a(w) represents the case 71 = 72 =1- It exhibits a minimum at uj =0 
and approaches unity for ui — > ±00. The corresponding correlation function in the left 
inset is negative for all times, i.e. only negative correlations are present (region ./V in 
Fig. la). By lowering 71 below one a maximum at uj — and two minima symmetric to 
uj =0 arise (dotted curve: 71 =0.4, 72 =1; full curve: 71 =0.2, 72 =1). The correlation 
function (full curve in the right inset of Fig. 2) belonging to the full curve a(uj) {a >1, 
P-region in Fig. la) has now two contributions with different signs: positive (dotted 
curve) and negative (dashed curve). The positive part of the correlation function is 
due to the bunching of tunneling events as discussed with respect to Fig. lb and the 
negative part is due to anti-bunching caused by Pauli's exclusion principle which is 
still present in the bunches of tunneling through the current carrying level. Therefore, 
the frequency-dependent Fano factor in Fig. 2 consists of a sum of two Lorentzians 
with either positive or negative sign and different FWHM corresponding to the two 
time scales related to the respective collector tunneling rates. 




Figure 2. Fano factor a vs. frequency w for 71 = 1 and 72 = 1 (dashed 
curve), 71 = 0.4 (dotted curve), and 71 = 0.2 (full curve). Insets: current-current 
correlation functions at the emitter barrier Cee vs. time t (full curves); dashed 
curve: negative part of Cee, dotted curve: positive part of Cee- 



4. Conclusions 

We discussed the effect of positive correlations in the tunneling current through two 
single-particle states in a regime where mutual Coulomb blocking leads to bunching 
of tunneling events. We showed that the current-current correlation function at one 
barrier consists of the sum of a negative and a positive term for an asymmetric choice 
of the tunneling rates of both QD states. Then, the corresponding spectral power 
density is a sum of two Lorentzians with different FWHM and different signs which 
gives rise to a maximum in the zero-frequency Fano factor and which can become 
larger than one (super-Poissonian shot noise). 
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